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; Abstract 

a 

^ ' Transport coefficients of gluon plasma are calculated for a SU(3) pure gauge model by lattice 

OC ■ QCD simulations on 16^ x 8 and 24^ x 8 lattices. Simulations are carried out at a slightly above the 

K^ . deconfinement transition temperature Tc^ where a new state of matter is currently being pursued 

^L \ in RHIC experiments. Our results show that the ratio of the shear viscosity to the entropy is less 

\^ ' than one and the bulk viscosity is consistent with zero in the region, 1.4 < T/Tc < 1.8 
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Introduction: RHIC experiments liave been bringing us many surprises. One of them is 
that the data are unexpectedly well described by the hydrodynamical model [l|. Experimen- 
tal data and phenomenological analyses suggest the quark-gluon plasma (QGP), or a new 
state of matter may be produced. See Ref. 12] for a review of RHIC experiments. Molnar 
and Gyulassy investigated the elliptic flow data using a Boltzmann-type equation for gluon 
scattering, and found that they needed a cross section about 50 times larger than expected 
in the perturbative QCD 3|. This indicates that the QGP state above the phase transition 
temperature, T^, is not a free gas of perturbed gluons. QCD-TARO collaboration measured 
the temporal meson propagators and found that their wave functions do not behave as free 
particles even at T ~ LST^; they conjectured that the strong interactions between the ther- 



mal gluons and quarks may provide binding forces. Recent 



more extended analyses of the 
and it was suggested that the 



temporal propagators were reported by three groups |^, |5, 
charmonium state survives until around 2Tc. 

The new state of matter produced at high temperatures in RHIC experiments is most 
likely not a weakly interacting plasma, but a strongly interacting quark-gluon system. In- 
vestigating the results in Ref. Isl], Teaney found that r]/ s ~ 0.04, where rf and s are the shear 
viscosity and the entropy, respectively 7]. Shuryak and Zahed have proposed a 'strongly 
coupled QGP' model for the new state of matter above TcSj, and argued that the QGP 
studied in RHIC is the most perfect fluid ever measured. Policastro et al. have calculated r] 
for the finite-temperature A/" = 4 supersymmetric Yang-Mills theory in the large A^, strong- 
coupling regime, and obtained t]/ s = l/47rJ2|. This value is found to be universal for theories 
with gravity duals and it is conjectured that rj/ s = l/An is a lower limit for all systems in 
nature [10|. 

It has now become highly desirable to study the nature of the quark gluon system, 
particulaly its hydrodynamical parameters such as the transport coefficients above T^ based 
on QCD in a non-perturbative manner. In this paper, we calculate the transport coefficients 
of QGP at a slightly above Tc, from the lattice simulations. Simulations are carried out in 
the quench approximation. For the calculation of the transport coefficients on a lattice, we 
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13l |. where the transport 



apply the formulation based on the linear response theory 
coefficients are calculated from Matsubara Green's function of energy momentum tensors. 
Numerical simulations of transport coefficients with this formulation were first carried out 
by Karsch and Wvld(l4l|. In their pioneering work, they performed the simulation on an 



8^ X 4 lattice but unfortunately, the size in the imaginary time direction was too small for 
the determination of the transport coefficients. 

We report here our simulation on a Nrp = 8 lattice with RG (renormalization group) 
improved action by Iwasaki. Our results are summarized as follows. 

1. The ratio of the shear viscosity to the entropy, r]/s, is small, i.e., less than one, but it 
is most probably larger than 1/Att. See Fig|21 

2. The bulk viscosity is less than the shear viscosity and is consistent with zero within 
the present statistics. 

3. For the heat conductivity, we could obtain no meaningful result. This is because, in 
pure gauge theory, there is no conserved current which transports the heat. J33| 

Preliminary results based on 16"^ x 8 and smaller lattices have been reported at lattice and 
Quark Matter conferences |l5l. Ilq . 

Transport CoefRcients in Linear Response Theory: The formulation for the transport 
coefficients of QGP in the framework of the linear response theory has been given in Refs. (ill . 
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13|. For the sake of consistency, we shall summarize the formula which will be used in 
the following calculations. 

Transport coefficients are calculated using the space-time integral of a retarded Green's 
function of energy momentum tensors, 

V=- J < Tuix, t)Ti2(f' , t') >ret, (1) 

\v + C = - J < Tnix,t)Tnix',t') >,,*, (2) 

X = -^ / < Toi(f,t)Toi(x',t') >ret, (3) 

where J = J d^x' J_ dtie^^^^^'^^ J_^ dt', and t], ( and x represent shear viscosity, bulk 
viscosity and heat conductivity, respectively. < T^^Tp^ >ret is the retarded Green's function 
of the energy momentum tensors at finite temperature. For the pure gauge theory, T^j^'s are 
written by the field strength tensors F^^,: 



Ff^i, are defined by plaquette variables on the lattice as U^^{x) = exp {ia^gF^u{x)). F^y are 
oMa.„ed eithe. by taHng the lo, of U,. d.ectly. o. by expand™ t/„„ «th .expect to a^,. 
In the following, we use the latter method to calculate -F^i/[1J|. |3J| 

It is difficult to calculate the retarded Green's functions in the lattice QCD, in which 
Matsubara Green's functions are measured. The retarded Green's functions are obtained by 
the analytic continuation. We obtain the numerical values of Matsubara Green's functions at 
discrete variables u;„ = 2TmT in the momentum space, while the retarded Green's functions 
are functions of the continuous variable po- Therefore, we need a bridge for the analytic 
continuation. 

Matsubara Green's functions Gp are expressed in a Fourier transformed form with the 
spectral function p: 

GMt) = Y.'--' h-^^^^- (5) 

^ J lUJn-OJ 

It is well known that the spectral function is common to both the retarded and Matsubara 

n 

Green's functions |17l|. The expression for the retarded Green's functions is obtained by 
putting u; — > po + i^- 

The determination of p(p, uj) is not straightforward, because in a numerical simulation, 
Matsubara Green's function has a finite number of points in the temperature direction, 
Nt/2. We must employ an ansatz for the spectral function with parameters, which are 
determined by fitting Matsubara Green's function. The simplest nontrivial ansatz for the 
spectral function has been proposed by Karsch and WyldjlJ], 

p{p = 0,uj) = -{- -— — ■ — -— — -), (6) 

TT (m — ujj^ + 7^ [171 + ujy + 7"^ 

where 7 represents the effects of interactions and is related to the imaginary part of the 



m- 



selfenergy. This ansatz is supported by perturbative calculations |11. 

Once we use this ansatz for the spectral function, the space time integral of the retarded 
Green's function can be calculated analytically. The result is 

27m ,„, 

where a represents the shear viscosity r^, bulk viscosity Q or heat conductivity x times T. 
At least three independent data points for Matsubara Green's functions are necessary to 
determine these parameters. 



In Ref. [ij], a simulation was carried out on a 8^ x 4 lattice, where two independent data 
points in the temperature direction are available. In this simulation, three parameters in 
the spectral function could not be determined. In order to determine A, 7 and m, we adopt 

Nt = S. 

Numerical Simulations: We calculate the transport coefficients in the SU{3) gauge the- 
ory for the regions a slightly above the transition temperature which are covered in RHIC 
experiments. We adopt Iwasaki's improved gauge action, which is closer to the renormal- 
ized trajectory than the plaquette action, and we obtain results close to the continuum limit 
on relatively coarse lattices 18|. We found that the fluctuation of the Matsubara Green's 
function is much suppressed comparing with the standard plaquette action 19| . 

We should ffist determine the critical (3 of Iwasaki's improved action on the Nt = 8 
lattice. For the Nr = 4 and 6 lattices, the critical (3 for this action were determined by the 






Tsukuba group |20|. We have carried out a simulation for (3c on a 16^ x 8 lattice; the results 
were reported in Ref.[l9|. However, the volume size was small, and we could obtain only a 
rough estimation of /3c, that is, 2.70 < (3c < 2.72. If we use the finite size scaling formula 
reported by the Tsukuba group, Pc at Nt = 8 becomes 2.72 < /3c < 2.74. The values of 
(3c determined by the simulation for A^^ = 4, 6, 8 do not yet satisfy the asymptotic two-loop 
scaling relation. We take (3 = 3.05, 3.2 and 3.3 as our simulation points. 

Matsubara Green's Function on N^ = 8 Lattice: The parameters of the simulations and the 
obtained statistics are summarized in Table 1. For Matsubara Green's functions Gn and 
G12 from which the shear and bulk viscosities are calculated, we can obtain reliable signals 
from approximately 0.8 x 10^ MC data on a 24'^ x 8 lattice. As an example, 6*12 is shown in 
FigII]for/3 = 3.3. 





(3 


total sweeps 


For equilibrium 


bin size 


163 


3.05 
3.2 
3.3 


1333900 
1212400 
1265500 


133900 
112100 
165500 


100000 
100000 
100000 


243 


3.05 
3.3 


861000 
784000 


61000 
84000 


100000 
100000 



TABLE I: Simulation parameters and statistics. Data at t = and t = 8 are not used for the fit. 



In the case of the 16'^ x 8 lattice, the errors are larger than the signal at r = 4, even with 
more than 10^ Monte Carlo (MC) data. The volume of 16^ may be too small for A^^ = 8. 

Gi4, from which the heat conductivity is calculated, has too large a background noise 
to extract a signal. Therefore, the fitting of Matsubara Green's function by the spectral 
function of Eq.ijH)) is carried out only for Gn and Gi2- 
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FIG. 1: Numerical data points and fitting results of Matsubara Green's function Gi2{t) on a 24^ x 8 
lattice 



Transport Coefficients of Gluon Plasma: The fitting of Matsubara Green's function by 
Eq.® was carried out by applying a non-linear least-square fitting program, SALS. Then 
the transport coefficients of the gluon plasma are calculated using Eq.Q. The errors are 
estimated by the jackknife method. After equilibrium is reached, the data are grouped into 
bins and the average of the data in each bin is treated as an independent data sample. The 
bin size is changed from 5 x 10^ to 1.2 x 10^. The results are independent of the bin size. 
In the following, the bin size is 100000 as shown in Table HI 





/3 


rja^ 


Ca' 


r?GeV3 


CGeV^ 


16^ 


3.05 
3.2 
3.3 


0.0018(28) 
0.0059(46) 
0.0043(90) 


-0.0015(29) 
-0.0025(20) 
-0.0041(142) 


0.054(82) 
0.281(223) 
0.283(590) 


-0.044(85) 
-0.122(90) 
-0.027(931) 


243 


3.05 
3.3 


0.0036(36) 
0.0072(30) 


-0.00095(288) 
-0.0031(26) 


0.106(108) 
0.471(194) 


-0.028(85) 
-0.201(167) 



TABLE II: Shear and bulk viscosities in non-dimension and in the physical units. The lattice 
scales, a~i = 3.09, 3.62 and 4.03 GeV for (3 = 3.05, 3.20 and 3.30, respectively. 

The results for the shear and bulk viscosities are given Table ini The bulk viscosity is equal 
to zero within error bars, while the shear viscosity remains finite. We do not see the size 
dependence. 

In the lattice calculations, the shear viscosity is calculated in the form rj x a^. In order 
to express it in physical units, we should know the lattice spacing a at each /3 value. For 
the estimation of a, we use the finite temperature transition point Pc- We take Pc = 2.73 for 
Nt = 8. The transition temperature is Tc = 276MeV\2Ql, and assume asymptotic two-loop 
scaling for the region /? > 2.73. The lattice spacing and the shear and bulk viscosities in 
the physical units are also listed in Table ITTl t]^^^ expressed in the physical units are slightly 
less than the ordinary hadron masses around T^. 

Entropy density: In a homogeneous system, the free energy has the form oi F = fV, and 

-/. Using the thermodynamical relation, U — TS 



then the pressure is p 
we obtain I21I 



s = S/V={e + p)/T, 



-TlogZ = F, 



(8) 



where e is the energy density. Using lattices with N^ = 8 and the integration method, 
p/T^fp^ = j^^ df3'N^{{S)T - {S)o), CP-PACS obtained p ande^ where {S)t and (^)o are the 
expectation values of the action density at temperature T[2^. We reconstruct the results 



from their numerical data of {S)t and {S)o and calculate the entropy density in Eq.( 

Concluding Remarks: In the high temperature lirnit, the transport coefficients have been 
calculated analytically by the perturbation method 11 
marized as follows. 




24 |2a 120, |22|. They are sum- 



(1) The bulk viscosity is smaller than the shear viscosity. This is consistent with our nu- 



merical results. (2) The shear viscosity in the next-to- leading- log is expressed by J27|, 
Vnll = (T^ / 9^)Ci/\og{fi* /mo) where niD = >/! + Nf/6gT, and for the pure gluon system 
Ci = 27.126 and fi*/T = 2.765 . 

There is a slight ambiguity in the relationship between coupling g and the temperature, 
and we use a simple form, g^'^ = 26olog(4T/A) with bo = llA^c/487r^. The scale parameter 
A on the lattice is set to be A/Tc ^1.5. For the entropy density, we use a hard-thermal loop 



result 



28l |. With these formulae, the perturbative rj/s can be compared with the results of 



numerical calculations. The result is shown in Fig|21 
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FIG. 2: The ratio of the shear viscosity to the entropy 



'Perturbative theory' is constructed from rj in Ref. 



as a function of T/Tc KSS bound is l/47r 



23] and s in Ref. 



2S|. 



Q. 



In this letter, we report the first lattice QCD result of the transport coefficients in the 
vicinity of the critical temperature. Although it still contains large errors, it may provide 
useful information for understanding QGP in these temperature regions. In particular, a 
small rj/s supports the success of the hydrodynamical description for QGP. Applicability 
conditions of the hydrodynamical model in quantum field theory were first considered in 
Ref. 29|. Together with experimental and phenomenological studies, the field theoretical 
approach will enrich our understanding of the new state of matter. We have shown here 
that the lattice QCD numerical simulations can provide useful information. 
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The next step is to obtain data with smaller systematic and statistical errors. If we can 
reduce the error bars in Fig|21 by a factor of two or three, we may realistically compare the 
data with the conjecture in Ref. [10[. We observed that Matsubara Green's function suffer 
from large fluctuations, but by using the improved action the fluctuations are significantly 
reduced. Another possibility for reducing the fluctuations may be to employ improved 
operators for T^,yJ30|. 

The results here depend on the ansatz of the spectral function of the Fourier transform 
of Matsubara Green's function. In order to test the functional form of the spectral function, 
we need more data points for Matsubara Green's function in the temperature direction, for 
which the most effective approach will be to apply an anisotropic lattice. If we have sufficient 
data points, the maximum entropy method is a promising way of determining the spectral 
function[31] which is free from the ansatz. Aarts and Resco pointed out, however, that it 
is difficult to extract transport coefficients in weakly- coupled theories from the euclidean 
lattice, since Green's function is insensitive to details of the spectral function p(uj) at small 

. New concepts will be necessary to overcome this diffculty. 
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